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In deterministic chaotic systems, the proba- 
bility that a trapped particle does not escape 
through some pre-specified leaking region or hole 
decays exponentially with time. One feature of 
this problem that has received recent attention 
is that even if the escape rate is of the order 
of the hole size it can exhibit strong fluctuations 
depending on the hole's position in the system's 
phase space. In this paper we compare the typ- 
ical escape rate obtained through two different 
procedures of introducing the holes in the setting 
of uniformly hyperbolic maps. In the first case, 
a small hole is placed at random in a specific lo- 
cation of the phase space and so the expected es- 
cape rate is the average over all possible hole po- 
sitions. In the second case, the random choice of 
hole position is performed independently at each 
time step and trajectory. This leads to a com- 
monly used estimate for the average escape rate 
and corresponds also to the physical picture that 
the map is partially leaking in the whole phase 
space, i.e. at every iterate of the map a fraction 
equal to the relative size of the hole is lost. While 
both these averages are equal in the limit of small 
holes, we show that for finite sized holes, the for- 
mer is greater than the latter and so escape is 
faster than expected. 



I. INTRODUCTION 

In recent years physical problems and mathematical 
results have motivated a renewed interest in the problem 
of placing holes through which trajectories can leak out 
from otherwise closed chaotic dynamical systems [1, 2] 
(for a recent review see Ref. [3]). The non-trivial aspect 
of this problem is that the properties of the open system 
depend sensitively on the position of the hole [4-12]. For 
instance, in Fig. 1 we show the escape rate 7j (i.e. the 
exponential rate of decay of smooth initial conditions) of 
the fully chaotic one-dimensional doubling map for dif- 



ferent positions i — . . . 2™ — 1 for holes of size h = 2~™ 
and n = 6. In the limit of small holes sizes h — > 0, it 
is well known that for any position (apart from a zero 
measure set) 7//1 — > 1. The most natural (and naive) 
approximation to finite (small) holes assumes that the 
(conditionally invariant) density remains uniform inside 
the open system and therefore 

j = -ln(l-h) = h+^h 2 + 0{h 3 ). (1) 

This expectation appears as horizontal lines in Fig. 1 
and correctly predicts the order of magnitude of the es- 
cape rate. The most striking deviation from this general 
feature are the deep minima, which are located at the 
positions of the lowest order periodic orbits [4]. Consid- 
ering a sequence of holes shrinking to a p-periodic point 
of the open doubling map, Ref. [8] shows rigorously that 
the escape rate is to first order given by 

^ = h{l-2- p )+o(h). (2) 

Observing this estimation in Fig. 1 we see that it provides 
an improvement over the naive estimation of Eq. (1) with 
7 < 7, and correctly captures some of the fluctuations in 
7, for finite size holes. 

In this paper we argue that there are many additional 
interesting features in the position dependence of 7,, 
apart from the minima, and that these problems require 
one to go beyond Eq. (2). To see this, consider the aver- 
age escape rate (7) = 2~™ Y^i=o 1 7« over the 2 n equally 
spaced hole positions (which correspond also to the ex- 
pected 7; if i is chosen randomly). From Eq. (2) we could 
expect that periodic orbits typically reduce the escape 
rate and therefore (7) would be smaller than 7. In this 
paper we will show that for the doubling map and for a 
large class of fully chaotic Markov systems the opposite 
is true and 

(7) > 7. (3) 

We will show that this inequality can be obtained through 
an asymptotic expansion to second order in hole size h 
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FIG. 1: (Color online) The escape rate 7 is plotted as a 
function of position for holes of size h — 2~ n and n — 6 
centered at ih + h/2 for i = 0, . . . , 2" — 1 . The true values are 
indicated by the black full curve, the first order approximation 
of Eq. (2) by the blue squares, while the red dots correspond 
to the second order approximations calculated in Eqs. (21) 
and (22). The width of each tower produced by the black full 
curve equals the size of the holes. The horizontal green line 
indicate the naive estimates of Eq. (1). 



of an exact periodic orbit formula (see Eq. (11) below). 
As seen in Fig. 1, this expansion provides a much better 
estimation of the true values of ji and can account for 
the many . 

Our paper is structured as follows: In Sec. II we intro- 
duce the formal definitions and background information 
needed for our calculations. In Sec. Ill we present the 
main results for the doubling map, including the new 
periodic orbit formula and the demonstration of the in- 
equality in Eq. (3). In Sec. IV we extend these results 
to other uniformly hyperbolic systems admitting a finite 
Markov partition including the two dimensional baker 
map. In Sec. V we discuss the case of skewed maps and 
show how the inequality can be generalized or even bro- 
ken. Finally, in Sec. VI we conclude and discuss the gen- 
erality of our results, possible extensions and applications 
to other research areas. 



II. FORMALITIES AND BACKGROUND 
MOTIVATION 

We consider the function p(x, n) describing the den- 
sity of representative points in phase space M. at time 
step n which evolves under the action of a deterministic 
closed map / : M. JA. The evolution of p(x, n) can be 
understood best in the language of operators such that 
p(x, n) — C n p(x, 0), where £ is the Perron- Frobenius 
operator associated with / [13]. Thus, the normalized 
invariant density p(x) of / is the eigenfunction related 
to the largest (in modulus) eigenvalue of C. In Ulam's 
method [14], C is approximated by an N x N trans- 



fer matrix T with elements corresponding to the transi- 
tion probabilities between the iV-partitioned phase space. 
There has been much work in recent years involving the 
convergence rates of Ulam's method [15] as N — > 00 as it 
can be used as a basis for rigorous computations [16]. For 
piecewise linear maps admitting a finite Markov partition 
such as the maps considered here, the leading eigenfunc- 
tion is piecewise constant, and so the treatment using a 
finite matrix is exact [13]. Thus, such maps are called 
Markov maps. 

A closed map can be opened by choosing the i th ele- 
ment of the partition as a hole Hi C M. through which 
trajectories can escape never to return. We denote / as 
the open map corresponding to /. At each iterate j of 
/, a proportion of mass Vj may be lost to Hi and there- 
fore / : M.\ Hi —> A4. Typically, in strongly chaotic 
maps (e.g. with exponential decay of correlations), the 
proportion of mass remaining in the system decays expo- 
nentially with n such that the escape rate 



7 = lim In P(n), 



(4) 



is well defined. Here, P(n) = n?=i(l — v j) * s cane d the 
survival probability. The leading eigenfunction of the 
open Perron-Frobenius operator C associated with / is 
now the conditionally invariant density p c (~x) (see also 



Ref. [17]) such that £> c (x) 



yn jO c (x), and the es- 



cape rate can be obtained from the leading eigenvalue 
A < 1 of C through the relation 7 = — In A (see also 
Refs. [18]). For Markov maps, p c (x) is again piecewise 
constant and so one looks at the leading eigenvalue of the 
finite matrix Tj = £, where the subscript i £ [0,N — 1] 
characterizes the position of Hi and is simply a zero- 
valued (i + l) th column. Note that the matrices T and 
Xi are non-negative and hence by the Perron-Frobenius 
theorem have a real positive eigenvalue (the Perron root) 
which is greater or equal in absolute value than all other 
eigenvalues. Moreover, they are sparse and hence one can 
utilize preconditioned iterative solvers for fast computa- 
tions. 

For the greater part of this paper we will use the 
one dimensional doubling map f(x) = 2x (mod 1) as 
our paradigm example and later generalize the results 
to other systems. We choose the doubling map as it is 
uniformly expanding and has invariant density p(x) = 1 
on the unit interval, thus making it amenable to analysis 
with the methods we employ. In particular, there is a 
close correspondence between the binary representation 
of a point x = 0. 01020,3 . . .2, and the symbolic dynamics 
for the partition {[0, 1/2), [1/2, 1]}, modulo minor details 
to do with dyadic rationals (i.e. fractions with denom- 
inator a power of 2), a zero measure set. Hence, with 

the partition I n = {I n ,i}i=o 1 wnere In,% = [M + l]2 _n 
for n > and i £ [0, 2" — 1], one can consider the open 
doubling map / with a Markov hole Hi = I„,i of size 
h = \Hi\ = 2~ n with respect to the relevant invariant 
measure of the closed map / (in this case Lebesgue). 
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The leading eigenfunction of the corresponding 2" x 2™ 
transfer matrix is piecewise constant on these intervals 
and so the leading eigenvalue can be calculated exactly 
as the root of a finite polynomial, or numerically to ar- 
bitrary precision. 

It originally came as a surprise when the escape rate of 
uniformly hyperbolic systems was shown to be strongly 
dependent on the position of Hi, allowing for the possibil- 
ity of escape through some holes to be as fast as through 
holes which are twice as big [4] . Similar results were also 
observed in two dimensional billiard models [19] (for a 
recent review see Ref. [20]). This striking observation 
(amongst others) was originally proved in Ref. [8] for a 
large class of hyperbolic maps, later and independently 
generalized in Ref. [6] , and was also extended and applied 
to the context of diffusion in Ref. [10]. Significantly, for 
the open doubling map Eq. (2) suggests that the local es- 
cape rate for each hole is lim ji/h = 1 — 2 _p , and since 

h— >o 

the number of aperiodic points (i.e. irrational numbers) 
is of full measure, almost every point in [0, 1] has a lo- 
cal escape rate equal to 1 [8]. Indeed, Eq. (2) seems to 
neglect the many local maxima for which Ji/h > 1 (see 
Fig. 1) and suggest that the 7 = h limit is approached 
from below since every finite sized hole contains a peri- 
odic orbit of finite period. 

A commonly used and very useful estimate for 7 is 
given by 7 as in Eq. (1), which is what Ref. [5] refers 
to as the "naive" estimate since it is equivalent to a bi- 
nomial estimate for recurrence times in a purely random 
process and also to the escape rate in the presence of 
strong noise [21]. Interestingly, a connection can also 
be made with random maps. We briefly recall that a 
random map is a discrete time process in which one of 
N maps is selected at random and applied separately at 
each time step [22]. Thus, in the spirit of Ref. [23], one 
can define an average transfer matrix for open Markov 

maps as T = X^o* li^ii with probability weights t^. 
The corresponding average escape rate is calculated from 
the leading eigenvalue of T. For the open doubling map 
(where N = 2 n ), if all but one qi are equal to zero, then 
we arc in the situation described previously with the es- 
cape rate being strongly position dependent. If however 
all transfer matrices are equally probable to occur then 
all q 4 = 2~ n , and T = (1- 2~ n )T where T is the transfer 
matrix for the closed system. This case corresponds to 
partial leakage (reflection coefficient equals to 1 — h) uni- 
form in the phase space and so the escape rate 7 equals 
that of Eq. (1) with h = 2~ n since T is measure preserv- 
ing and so has largest eigenvalue equal to 1. 

We are interested in the following natural question: 
If a Markov hole is chosen at random, what escape rate 
should one expect? or equivalently, What is the average 
of the curve shown in Fig. 1 ? We thus define this average 
for a Markov map with N possible hole positions as 



N-i 



For the doubling map, all hi = 2~" and therefore (7) is 
just an arithmetic mean. In the limit of large n we have 
that (7) converges to the size of the hole in agreement 
with both Eqs. (1) and (2). It is the intention of this 
paper however to show that for a large class of Markov 
maps this limit is actually approached from above rather 
than below. We do this by first deriving an exact peri- 
odic orbit formula for Markov maps (Sec. Ill A) which 
we then asymptotically expand to second order in hole 
size (Sec. IIIB). We then take the arithmetic mean of all 
2 n escape rates and obtain an asymptotic expansion of 
(7) (Sec. IIIC) leading to the inequality (3) for n > 1. 
Thus the escape rate is typically faster than the expected 
naive estimate 7. 



III. MAIN RESULTS 

A. Periodic orbit formula for Markov maps 

As discussed above, the leading eigenvalue A for the 
open doubling map is given by the solution of a polyno- 
mial equation, the characteristic equation of the transfer 
matrix Tj . To make progress we need an explicit expres- 
sion for this polynomial. While this has been discussed 
in a number of contexts, the form in which we express 
this, a new periodic orbit formula, is of interest in its own 
right. 

Historically, this problem has been considered from the 
point of view of the waiting time distribution for finding 
a given fixed sequence S of n symbols in a sequence built 
drawing an independent identically distributed random 
variable, which has practical applications in computer 
search algorithms and DNA sequence analysis. While 
similar expressions were stated as early as 1966 [24], the 
most convenient starting point is Thm 2.1 of Ref. [25] 
which gives (using our notation) the following recursion 
relation for the probability Wt of stopping after exactly 
t > n symbols: 



w t =2- n -2- n J2^ P - E w v 2-^ Xp +n-u 



p—n p— t— n+1 



(6) 



where \r is one if the first and last r symbols of S are 
identical, otherwise zero. 

We expect w t to decay exponentially for large t such 
that wt = cA*(l+o(i~ Q )) for any power a > 0. Moreover, 
since the wt are probabilities they sum to one. Thus, by 
the use of the formula for a geometric series, the first two 
terms on the RHS become 



2- 1-£>P =2- £ 



p—t— n+1 



<7> = E 



(5) 



(2A)" 



(7) 



2=0 



1-A 



-cA-*(l + o(t- Q )). 
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Substituting this into Eq. (6) , writing J = t — p in last 
sum on the RHS of Eq. (6) and relabeling J as p, dividing 
everything by cA* and then taking the limit t — > oo we 
arrive at an exact formula 



A 



1 - A 



(2A)" 



n-l 

E 

P =i 



(2X)- P Xr, 



(8) 



We now connect this problem to the characteristic 
equation of the transfer matrix Ti (not restricted to the 
Ti of the doubling map). First we notice that for each 
hole indexed by i = . . . 2 n — 1, there is a fixed point of 
the map f n . We then interpret the symbolic sequence 
Si of the periodic orbit associated to this fixed point as 
the symbolic sequence of the hole Hi. Significantly, in 
the case of the doubling map, Si is precisely the n-digit 
binary representation of the hole index i. Next we no- 
tice that Xn-p used above indicates exactly the number 
(zero or one) of periodic points of length p in the interval 
Hi = I 1l: i. Thus Eq. (8) can be written as 



1 - A 



(2A)"" 



n-l 



E 

P=l x:/ p (x)=x 



(2A)- p - 



(9) 



A more compact version of this expression can be ob- 
tained by noting that each period p > n has exactly 2 P ~" 
periodic points in the interval Hi irrespective of the value 
of i, and allow for a divergent sum, interpreted using the 
formula for a geometric series such that 

oo OO . 

J2 E (2A)-* = X>- n (2A)-' = _ I (2A) -n 

P="x:/ P (x)=x P=n 



(10) 



Therefore, we finally arrive at our first main result 



! + E E (2A)- p = 0. 



p=l x:/ p (x)=x 



Note that this differs from other periodic expansions in 
the literature [26] in that it enumerates periodic points 
in the hole rather than periodic orbits that avoid it (see 
also [5]), the sum is over all periodic points at all periods, 
and it is always divergent. The derivation naturally gen- 
eralizes to Markov maps with more than two branches, 
including where each branch has a different expansion 
factor (i.e. skewed maps); here the factor 2 _p is replaced 
by the expansion factor of the relevant periodic orbit. 



B. Small hole asymptotics 



For the doubling map with Markov holes Hi , all of size 
h = 2~ n , Eq. (11) is exactly equal to the characteristic 
equation of Tj multiplied by hX/(X— 1). While the sum in 
Eq. (11) diverges, it can be interpreted as a geometric se- 
ries and hence Eq. (9) can be expressed in the convergent 
form 



(A - 1) ( 2™A n ~ 1 + 



-p A "-p-1) +i = o, (12) 



where V C [1, n—l] is the set of periods p < n for periodic 
points in Hi. 

We now expand the escape rate in powers of the hole 
size h; as can be seen from a similar calculation [27] this 
can be useful whether or not the function is smooth. We 
write 7 = + 7^' ft, 2 + . . ., so that the small hole 

expansion for the corresponding eigenvalue is of the form 
A = 1 — A^fr, — \^h 2 up to second order such that 
7 « = A« and 7 < 2 ) = (A^ 2 ) 4- (A«) 2 /2). Substituting 
(11) these into Eq. (12) we obtain 



E 

pev 



(l - („ - p - l)A«ft -((n-p- 1)A( 2 ) - (n-p-l)(n-p-2) {x{1))2 \ ^ _ \ 



(13) 



4-1 = 0, 



where Vq = V U {0}. Collecting terms of order h° = 1 and h 1 = 2 ™ gives 



„<4) 



E^ 

vPS-Po 



(14) 
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(a) _ Ep 67 ,„("-P- 1/2)2-' 

' / \ 3 ' 



(j2 P ev 



2-p 



(15) 



respectively. The expansion can be performed nicely to 
arbitrary order. 

We now make the important observation, that each 
Markov hole may be one of two types. Type A holes con- 
tain a single primitive periodic orbit of period p <E [1, n/2] 
which is repeated up to n — 1, and may also contain pe- 
riodic orbits with periods (n/2, n — 1] which are non- 
repeats of p. Typically, type A holes are associated with 
escape rates 7, < h due to the short periodic orbits. Type 
B holes contain none or many primitive periodic orbits 
of periods p £ (n/2, n — 1] and usually have escape rates 
li ?i h. In particular, the maximal escape rates appear 
in type B holes and the number of such holes is given 
by the integer sequence [28] A003000 which refers to the 
number of "bifix-free" words of length n over a two-letter 
alphabet and increases like ~ 0.2678 x 2™. We call the 
sets of type A and B holes A, B C [1,2™] respectively 
such that |^4| + \B\ = 2™. This simple classification of 
holes turns out to efficiently capture the fluctuations of 
7i and is therefore key in calculating the average given 
by Eq. (5). We will now treat each case separately. 



have that 



rc -1/2 p2-f 
1 - 2-p ~ (1 - 2-p) 2 



p:p\p 
p£V 

E 

p:p\p 

per 



" ~ l ' ~ 2 ' 2 



" "' '" ~ 2 ' 2 



(18) 



Approximating the denominator of 7^ by (1 — 2- p )- 3 , 
ignoring the sum (which is of higher order) and combin- 
ing with (18), we obtain that 



7< 2 > = ( n - ~ ) (1 - 2- p ) 2 - p2~ p (l - 2- fJ '« 



2. Type B holes 



(19) 



As there may be many primitive periodic orbits of 
length p > n/2 we expand Eq. (14) in a binomial series 
to get 



Type A holes 



As there is only a single primitive periodic orbit of 
period p < n/2, the sums in (14) and (15) must count 
repeats of this orbit up to n — 1, that is 



2~ p = J2 2 ~ tp + E 2 ~ p 

pev i=i p-p\p 

pev 

1 _ 9— ,n P 

2~ p - ^— + 2 ~ p > 

p:p\p 

pev 



(16) 



1 - 2" 



where p = min('P) < n/2 and m — L^^rJ — 1- Substi- 
tuting back into (14) we obtain that 



p:p\p 

pev 



(17) 



Notice that the leading order term is equal to Eq. (2) 
since the sum in (17) is at most of order ~ 2"™ /2 . 

We now consider the numerator of 7W in (15). We 



v(!) 



i + E 2 ~ p 



pev 




(20) 



For 7^ 2 ) it is sufficient to keep only the leading order term 
given by p = such that 7^' = n — 1/2. 



3. Comparison with true escape rate 

We define 74 as the escape rate corresponding to hole 
Hi = I n ,i such that we have the following approximations 
for type A holes 

lieA = (l - 2- p + 2-(™ +1 )p - 2 ~ P ) h 

p:p\p 

pev 

(n - - 2-P) 2 - p2"P(l - 2-P)] e + o(h 2 ), 



and for type B holes 



lieB 



l-£2- 

pev 



1 

n 

2 



(21) 



/i 2 + o(/i 2 ). (22) 
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Fig. 1 shows a comparison of Eq. (2) and Eqs. (21) and 
(22) against the true values of 7 for n = 6, 7. It is clear 
that the improved asymptotics now successfully capture 
both local maxima and minima of the escape rate. 



C. The average escape rate 

Equipped with Eqs. (21) and (22) we can now take 
the arithmetic mean of all j^s to obtain an asymptotic 
approximation to (7). To do this we need to know the 
number of periodic points a(p) in [0, 1] with primitive 
orbit lengths p > 1. For the doubling map / , we have 
that a(p) = 2, 2, 6, 12, 30, 54, . . . which exactly equals the 
number of aperiodic binary strings of length p and is 
given by the integer sequence ^4027375 described by the 
formula[28] 



p/d 



d\p 



(23) 



where /i(<i) here is the number theoretic Mobius function. 
Hence, there are \A\ = Ylp=i a (p) tyP e A holes and \B\ = 
2" - |.4| type B holes. 

We first sum over the escape rates of type A holes. For 
this it is sufficient to keep just the leading order term of 
jieA = (1 - 2~ p )h. We have that 



1/2 



ieA 



It = 



o(p)2" 



n/2 n/(2d) 



™/ 2 2 - fc (' i - 1 ) - 



d=2 



2 -L»/4j + 



2- 2 L«/6j 



d=2 




(24) 



where we have set p = jd and k — [n/(2d)\ in order to 
exchange the order of the sums in the second equality and 
[x\ is the integer part of x (floor function). Therefore, 
(24) converges exponentially to 



li = [\A\- 



ieA 
n/2 - M (rf) 



(25) 



1.382714 for n > 1. 



where K = J2 d =2 i-2<*-i 

We now sum over the escape rates of type B holes. 
Here it is necessary to keep terms up to second order in 




FIG. 2: (Color online) The difference between (7} and 7 was 
computed numerically using long double (19 digit) precision 
for n < 20 (blue dots). The red dashed line is given by the 
prediction (7) — 7 = «2~ 2n . 



h. We have that 



£-* = E 



\h \ \n--\ir 



B\ 



= {\B\- 



\ pev J 
E a^y+in-^B^ (26) 
-l) - [2-L«/*J _.. V\ h 



p=L»i/2j+l 

n 



Therefore, (26) converges exponentially to 

I>,= (isi-[ nl 1 



(27) 



since \B\ - 2™ - 2f +1 for n > 1. 

Using (25) and (27) we may now calculate the average 
(mean) escape rate 



2™-l 



<7>=2- n 5>i = 2- tt X^ + E-fl 

i=0 \ieA ieB 

V 



(28) 



, +■ - , 2- 



Since k > 0, the above calculations show that for small 
but finite Markov holes, 7 < (7) with equality only in 
the limit of n — > 00. The difference between the two is 
(7) — 7 = k2~ 2 " and is shown in Fig. 2 on a log-linear 
scale. Even if the derivation of Eq. (28) does not provide 
a rigorous proof of inequality (3) - we have not prop- 
erly bound Eqs. (21) and (22) - the excellent numerical 
agreement for n > 3 observed in Fig. 2 strongly suggests 
that fluctuation due to higher order terms are negligible. 

We now briefly summarize the results presented in this 
section. We first obtained a periodic orbit formula (11) 
for the leading eigenvalue A of the open transfer matrix Ti 
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FIG. 3: (Color online) Plot of jk for the baker map (30) for Markov holes h — \I m , n ,i,j\ for m — n — 1 and n = 2, 3, 4 and 5. 
The colors have been scaled in each case so that the minimum and maximum escape rates are represented by black and white 
boxes respectively (see color scale on the right). 



of uniformly hyperbolic maps admitting a finite Markov 
partition. Using the doubling map as our main example 
we have divided all holes into two distinct groups which 
correspond to escape rates typically greater or smaller 
than the expected naive estimate (1), and also obtained 
improved asymptotic formulas for them (see Eqs. (21) 
and (22)) to second order in hole size. Taking the arith- 
metic mean of all different escape rates, we attained an 
asymptotic expansion for (7), thus analytically showing 
that escape in the binary shift is faster than expected (see 
Eq. (3)) for finite size holes. We now consider generaliza- 
tions of these results to other open dynamical systems. 



IV. GENERALIZATIONS AND EXTENSIONS 
TO OTHER MAPS 

A. Linear expanding maps 

The above results can be easily generalized for linear 
expanding maps on the interval with uniform invariant 
densities of the form x 1— > sx (mod 1), s£Z and \s\ > 1. 
This can be done by using \s\ instead of 2 in Eq. (12), and 
expanding A in terms of the size of the new Markov holes 
h = \s\~ n . We remark that maps with s < have been 
recently used in the context of polygonal billiards with 
(non-conservative) "pinball" type dynamics to rigorously 
prove hyperbolicity [29]. More specifically, the case of 
s = — 2 corresponds to the so called "slap" map of an 
equilateral triangle billiard. 



B. Tent map 



The tent map 

fix) = 1-2 



, for x G [0, 1] 



(29) 



stretches [0, 1] to twice its original length and then folds 
it in half back onto [0, 1]. Although the map has a uni- 
form invariant density, the dynamics no longer commutes 



with the symmetry x — > 1 — x. Nevertheless, our results 
also apply here since the tent map is a metric conjugacy 
onto the left shift symbolic space and so shares the same 
binary symbolic dynamics and hierarchy of periodic se- 
quences as the doubling map. 



C. Baker map 



The two dimensional baker map 

((2x,y/2), for < x < 1/2, 

I{X ' y > |(2-2x, l-y/2), for 1/2<X<1, 



(30) 



is area preserving with respect to Lebesgue (i.e. p(x, y) = 
1), and is a two-dimensional analog of the tent map. The 
unit square is squeezed uniformly two times in the ver- 
tical y direction and stretched in the horizontal x direc- 
tion. It is then cut in half, and the right half folded over 
and placed on top of the left half. Hence, it is topo- 
logically conjugate to the Smale horseshoe map. Unlike 
all previously mentioned maps however the baker map is 
invertible and so can mimic chaotic dynamics in Hamil- 
tonian systems. The quantum version of the map has 
been used to explore the classical to quantum correspon- 
dence in the semiclassical limit [30], and to study the 
emergence of fractal Weyl laws in the quantum theory of 
open systems [31]. 

As with linearly expanding maps, the binary repre- 
sentations of points x and y in [0,1] 2 have a close cor- 
respondence with the symbolic dynamics for the parti- 
tion {[0, 1/2), [1/2, 1]} x {[0, 1/2), [1/2, 1]}. Therefore the 
Baker map also shares the same hierarchy of periodic se- 
quences as the doubling map. Moreover, the eigenfunc- 
tions of C n are piecewise constant on the rectangles given 
by 



i + 1 



2™ ' 2" 



(31) 



for i € [0,2 m - 1], j £ [0,2™ - 1], and n,m e N+. The 
map is Markov and thus the evolution of densities can be 
exactly described by a 2 n+m x 2 n+m transfer matrix T. 



The baker map is opened by choosing a Markov hole 
Hi j = I m ,na,j of size h = 2 _ (' i+m ) and setting the k th 
column of T equal to zero and denoting the open transfer 
matrix as T^. Fig. 3 shows the variation of 7^ as a 
function of position for m — n — 1 and n = 2,3, 4, and 
5. Since the leading eigenvalue of Tj. is a solution to 
Eq. (11), the small hole asymptotic formulas (21) and 
(22) of section 4.3 hold. Moreover, for n 3> 1 we have that 
(7)— 7 = n2~ 2 ( m+n \ We have confirmed this numerically 
for n < 10. 



D. Logistic map 



0.0 0.2 



There also exist examples with nonlinear dynamical 
equations and non-uniform invariant densities for which 
the above results hold. One such example is the logistic 
map 

f(x) = 4x(l -x), for x e [0, 1] (32) 



FIG. 4: Escape rate as a function of hole position for Markov 
holes Hi of different sizes hi for n = 4 (blue) and n — 5 (red, 
dashed) for the skewed map Eq. (34) with u — 3. The holes 
are centered at the dots and are of size equal to the corre- 
sponding tower widths. The horizontal (black) lines corre- 
spond to 71. 



with a partition of the interval given by 

2 f'(*+l)w 



sm 



2 n+l 



, sm 



2 n+l 



(33) 



for i — 0...2" — 1, and holes Hi — I ni . This is 
because the logistic map and the tent map are met- 
rically conjugate through the nonlinear transformation 
y = sin 2 (ttx/2). The average escape rate is given by 
Eq. (5) with N = 2 n and hi = 2~ n since the invariant 
density of the (closed) logistic map is p(x) = — , 1 

ttW :c(J x) 



V. EXPANDING MAPS WITH DISTORTION 

Here we argue that a sufficient requirement for inequal- 
ity (3) to hold is that the dynamics of / or some metric 
conjugacy of / are uniformly hyperbolic with a constant 
expansion rate. As a simple example of a case with vari- 
able expansion rate and the opposite inequality, consider 



- [ ux, < x < l/u 

)v(x — l/u), l/u<x<l, 



(34) 



where v = u/{u — 1), u > 1, and u 7^ 2. Map (34) 
has a uniform invariant density p(x) = 1 on the unit 
interval and shares the same binary symbolic dynamics 
as the doubling map for the partition {[0, l/u), [l/u, 1]}, 
and hence hierarchy of periodic sequences. However, the 
map does not have a constant expansion rate and its 
Lyapunov exponent is given by ^ In it 4- ^ In?; . 

The natural partition of [0, 1] under the dynamics of 
/ is composed of intervals of variable lengths \I n ,i\ € 



-(n-l),,-l 



I £ [0,n]} with binomial occurrences. Hence 



the corresponding Markov holes Hi are not all of the same 
size. The escape rates for the case of u = 3 are plotted in 



Fig. 4 for n = 4,5. Notice that there is only one minimal 
escape rate at Hq. 

The generalization of Eq. (2) for skewed maps (includ- 
ing skewed tent maps) is given by [6] 



j = h (l-A~ l ) +o(hi), 



(35) 



where A x is the stability eigenvalue of the periodic or- 
bit starting from x € Hi, and is infinite if the orbit of 
x is aperiodic. Hence, the stability of a periodic orbit 
depends not only on its period but also on its symbolic 
sequence. This implies that while the periodic orbit for- 
mula of Eq. (11) can be naturally generalized to the cur- 
rent setting by replacing 2 P by A x , the results of Sec IIIB 
need to be reformulated in terms of stability orderings 
rather than just the period lengths. 

For open skewed maps, one also needs to reconsider 
the "naive" estimate (1), as there now may exist sev- 
eral possible candidates. We concentrate on the fol- 
lowing three: (i) Clearly, the average hole size is 2~™ 
(independently of u) and so a first candidate is simply 

71 = — ln(l — 2~"). (ii) Assuming that the condition- 
ally invariant density remains constant inside the open 
system, we obtain a second estimate which we denote as 

72 = - Ei=o 1 fh ln(l - In) = u- 2n {u - 2u + 2)™ + . . . for 
n ^> 1, where we have used the fact that there are (?) 
holes of size u~( n ~^v~ l . (iii) A third and not so naive es- 
timate for the expected escape rate can be obtained from 
the random map setting, discussed in Sec. II, by looking 
at the leading eigenvalue of T = J2i=o hiTi. For the 
case of n = 2 this can be expressed in closed form as 



3(u-l) + ^(u~l)(Au 2 -7u+T) \ 
73 = -ln — 2 . (36) 
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FIG. 5: (Color online) A comparison between the average 
escape rate (7} and the three naive estimates (71,72,73) as a 
function of u. The difference between the averages is plotted 
for n — 2,3, and 4, corresponding to the full, dashed and 
dotted lines respectively. Notice that for 72 and 73, inequality 
(3) can be broken. 



Note that estimates 71 , 72 and 73 are smooth functions of 
u and for u = 2 are equal to the original naive estimate 
(1). Furthermore, we remark that the three estimates 
described above, contain increasingly more information 
about the system at hand with 71 considering only the 
number of possible holes, 72 also considering the variation 
in hole sizes, and finally 73 also capturing some of the 
dynamics through the average transfer matrix T. 

A comparison between the exact value of (7) obtained 
numerically from (5) and the three candidate "naive" es- 
timates is shown in Fig. 5 as a function of u for n = 2, 3, 
and 4. From the figure, it is clearly seen that the inequal- 
ity is satisfied for case (i) entailing 71 < (7) but not for 
cases (ii) and (iii) such that 72,3 > (7) for sufficient dis- 
tortion, i.e. when u/v > 1 or m/d < 1. fn contrast, the 
inequality in case (i) is strengthened with added distor- 
tion. Therefore, the benchmark for comparison with (7) 
may be of equal importance as the underlying dynamics 
of the map considered; it might be interesting to investi- 
gate these further. An alternative naive estimate would 
be to consider holes of constant size, i.e. non-Markov 
holes. 

We briefly discuss the kinks observed in Fig. 5. For suf- 
ficiently large it, the eigenvalue associated with the fixed 
point at 1 (equal to 1/v), increases until it is equal to 
the eigenvalue associated with the whole repeller in the 
rest of the dynamics. This degeneracy can only occur 
for holes which can break the transitivity of the dynam- 
ics and hence cause competing escape rates, in this case 
crossing over at the parameter value u — n. By symme- 
try, the same occurs at u = (i- e v = n). A similar 



effect was observed for the hyperbolic stadium billiard 
[9] with two holes corresponding to a source and a sink 
which can be placed in such a way as to produce asym- 
metric transport. In both cases, the presence of the hole 
renders the dynamics non-transitive. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper we studied the dependence of the escape 
rate in uniformly hyperbolic dynamical systems on the 
location of Markov holes, with focus on the "typical" es- 
cape rate. We have analytically derived an exact periodic 
orbit formula which provides the escape rate as a function 
of the periodic orbits inside the hole. Using an asymp- 
totic expansion to second order in the hole size allowed 
us to obtain an expression for the average escape rate 
(7) which is larger than the expectation 7 from the size 
of the hole alone, a surprising result in view of the fact 
that short periodic orbits typically lead to 7 < 7. These 
results were illustrated and complemented with an exact 
numerical analysis and are valid for systems conjugate to 
the binary shift (e.g. linearly expanding maps, tent map, 
logistic map, and baker map) 

An important question is the generality of these results 
to different classes of systems. In this regard, the skewed 
map discussed in Sec. V shows that when the Markov 
partitions have different sizes, there are different possible 
generalizations of the naive estimate 7 and that our main 
inequality (7) > 7 holds only in one of the cases. We 
have also performed numerical simulations in the chaotic 
diamond billiard, which is beyond the class of systems 
investigated here as there is no finite Markov partition of 
the phase space. Our observations (not shown) indicate 
that if the hole is placed along the border of the billiard 
the inverse inequality i.e. (7) < 7 appears to hold for all 
hole sizes. In general, it is an interesting open problem to 
verify in which classes of systems the inequality (or the 
reversed inequality) holds systematically for all (small) 
holes sizes. 

Another assumption in our analysis which asks for gen- 
eralization is the choice of the holes to coincide with the 
Markov partitions of the maps. It is natural to consider 
also the doubling map with a non-Markov hole of size 
h = 2~" centered at any point Xh € [h/2, 1 — h/2]. When 
Xh is varied smoothly, the escape rate j Xh becomes a 
highly non-smooth function which in the limit of h — > 
contains a dense set of maxima and minima; a fractal 
function [27]. For finite hole size however, j Xh has many 
locally constant intervals (see also Ref. [II]) which are 
well described by the discrete 7j's corresponding to the 2 n 
Markov holes, modulo fluctuations of order h 2 \nh [27]. 
This happens because in the doubling map the largest 
invariant set that never reaches the hole is a locally con- 
stant function of the hole size, shape and position [II]. In 
Fig. 6 we have confirmed this picture by varying contin- 
uously the position of the hole of size h — 2™ (we expect 
similar results also for size h =^ 2"). Our basic inequality 
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FIG. 6: (Color online) Top: The escape rate j Xh for a hole 
of size h — 2~ n (for n = 4,5) centered at ij, € [ft/2, 1 — 
ft/2] (red jagged line) is compared with the escape rate 7; 
through 2 n Markov holes (blue line connecting dots). The 
horizontal distance between two dots is equal to h. Bottom: 
Plot of \j Xh — 7»| for a hole of size ft = 2 -5 . The discrete 
values of 7; seem to capture the main character of r y Xh , modulo 
fluctuations of order ft 2 In ft. 



(3) remains valid if the average (7) is computed using 
the continuum of possible hole positions (instead of only 
those at the Markov partitions). Altogether, these results 



suggest that the Markov holes forms a "Markov skeleton" 
of the full curve and can be considered representative of 
the more elaborate problem. 

Finally, it is worth comparing our results to different 
approaches. In Ref. [10] it was shown analytically that 
the average diffusion coefficient, calculated over different 
hole positions in a related dynamical system, is exactly 
equal to the size of the holes h. In contrast, our re- 
sults show that the escape rate is always larger than h. 
Ref. [32] made the insightful suggestion to interpret the 
transfer matrix T of our Markov maps as the adjacency 
matrix (also called the connection matrix) of a network 
(or weighted directional graph). In this analogy, making 
a hole in the phase space M. is equivalent to consider 
one of the N network vertices to be a "sink" (not able 
to transmit information) [7, 32, 33]. Our results show 
that, for the strongly connected networks corresponding 
to Markov maps, the random choice of a sink leads to 
a loss of information which decays typically faster than 
the expected e~ x l N . It is also worth noting that, beyond 
the interpretation mentioned in the first paragraph of this 
paper, there are many different alternative approaches to 
random leaks such as the randomly perturbed metastable 
interval maps considered in Ref. [34]. 
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